A periodic elastic medium in which periodicity is relevant 
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We analyze, in both (1 + 1)- and (2 + 1)- dimensions, a periodic elastic medium in which the 
periodicity is such that at long distances the behavior is always in the random-substrate universality 
class. This contrasts with the models with an additive periodic potential in which, according to the 
field theoretic analysis of Bouchaud and Georges and more recently of Emig and Nattermann, 
the random manifold class dominates at long distances in (1 + 1)- and (2 + l)-dimensions. The 
models we use are random-bond Ising interfaces in hypercubic lattices. The exchange constants 
are random in a slab of size L"^"^ x A and these coupling constants are periodically repeated along 
either {10} or {11} (in (1 -I- l)-dimensions) and {100} or {111} (in (2 -I- l)-dimensions). Exact 
ground-state calculations confirm scaling arguments which predict that the surface roughness w 
behaves as: w ~ L'^^^, L <^ Lc and w ~ L^''^, L ^ Lc, with Lc ^ X^''^ in (1 -I- l)-dimensions and; 
w ~ L° *^, L < Lc and w ~ ln(L), L > Lc, with Lc ^ A^'^® in (2 -I- l)-dimensions. 
PACS # 05.70.Np, 75.10.Nr, 02.60.Pn, 68.35.01 



I. INTRODUCTION 

Periodic elastic media arise in a surprising array of 
problems, including spin or charge density waves; flux 
line lattices; and random magnets. A model frequently 
used to describe a manifold, defined by the single 

valued height variable in a periodic elastic medium 
is 

Hpem = I df[^[Vhir)]' + rj[h{f)] + V,[hir)]} , (1) 

where is a periodic potential in the height direction 
and the random potential rj is not periodic. This is di- 
rectly analogous to the model used to study lattice ef- 
fects in thermal roughening and in field theoretic studies 
of commensurate phases in Ising magnets with compet- 
ing interactions. In the model (p, the periodic potential 
is non-random and tends to pin the interface while the 
quenched random pinning ri[h(r)] prefers to make the 
interface wander. The surface tension term ^[V/i(r)]^ 
seeks a flat interface and also competes with the quenched 
random pinning. Field theoretic calculations [0-^ sug- 
gest that at long distances, for (1-1-1)- and (2-1-1)- dimen- 
sional interfaces, the periodic pinning potential is irrel- 
evant, and hence that the interface scaling behavior is 
in the random-bond-Ising universality class where width 
= (ft,^) — (h)^ ~ L^'' with the roughness exponent 
C = 2/3 in (1 + 1), and C ~ 0.21(4 - d) in (d+l), d>2 
|0-p| . Note that lattice calculations are strongly affected 
by a lattice pinning potential and have a flat phase even 
for large lattice sizes [||. 

Another problem which has been heavily studied is the 
random substrate problem |£-11|. This was introduced 
to model the effect of a random substrate on layers of 
absorbed atoms, and also serves as a model for the ef- 
fect of a p-fold random field on the XY-model M . There 



is now a consensus that there is a disorder dominated 
glassy phase in this model (in 2 substrate dimensions) at 
low temperatures that is reflected in long distance cor- 
relations which behave as C(r) ~ In^ |r| (in contrast to 
thermally rough correlations in dimension (2-1- 1), which 
grow as as C(r) ~ lii|r|). There has been some uncer- 
tainty about whether the leading order correlations found 
by Cardy and Ostlund (CO) |^ are correct, with func- 
tional renormalization group calculations agreeing with 
CO 1 10 1^, and variational calculations disagreeing. The 
substrate roughness is randomly drawn from the interval 
(0,1) (in lattice units). This corresponds to a different 
sort of periodic elastic medium than that described in 
(|l|) above. Here, the random substrate leads to a period- 
ically repeated disorder seen by an interface lying above 
the random substrate. This arises due to the fact that 
the first, third, fifth, etc. atoms deposited at the same 
position on the random substrate see exactly the same 
disorder when they land. This corresponds to a random- 
bond Ising magnet in which the disorder is repeated with 
period A = 2 along the growth direction. In general, the 
disorder may range over a scale (0, A — 1), and this leads 
to a periodic variation in the disorder on length scale A. 
The continuum model for this system is simply. 



df 



{l[Vh{r)r + rj[h{f)]}, (2) 



but where rj is periodic in h{r), so that we require 
rj[h[f) -I- A] = rj[h{r)\. There has been considerable study 
of the random substrate (A = 2) problem, with the early 
controversy now being resolved in favor of a "super- 
rough" "Bragg-glass" phase in (2 -I- l)-dimensions in 
which w ^ ln(L). Exact ground-state calculations have 
been very useful in resolving this controversy [^-|l5|. It 
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is quite easy to see (see Section III) that in (1 + 1)- 
dimensions, the random substrate problem behaves as 
a random walk (RW), so that w L^/^. Note however 
it has been recently argued that although typical dislo- 
cations do not destroy the "Bragg-glass" ground state, 
optimal dislocations have negative energy, and hence are 
expected to destroy the Bragg glass in (2 + l)-dimcnsions 

In this paper we study the Hamiltonian (||) as a func- 
tion of the periodicity A of the disorder. We show that 
at long length scales in (1 -I- 1)- and (2 + l)-dimensions, 
the periodicity is relevant and that the random substrate 
universality class holds. The paper is arranged as fol- 
lows: Section || sets up the model and describes the way 
in which we calculate the exact positions of interfaces 
in random Ising magnets. The scaling theory describing 
the behavior of these interfaces is developed and tested 
in Section III . We give a brief conclusion in Section |v| . 



II. 



DISCRETE MODEL AND EXACT 
ALGORITHM 



The model which we use to analyze the effect of peri- 
odic disorder on interface properties is a spin-half Ising 
system with random bonds (RB) on square and cubic 
lattices. The Hamiltonian is given by. 



"Hrb — — JijSiSj 



(3) 



where > are coupling constants and the spin vari- 
ables Si take the values ±1. The spins on two opposite 
boundaries of the lattices, z = 1 and z = L, are fixed 
and have opposite signs so that an interface must exist 
in the lattice. Our calculations are at zero temperature 
and we find the ground-state interface properties for in- 
terfaces whose average normals lie in the {10} or {11} 
directions of square lattices and in the {100} or {111} 
directions of cubic lattices. The coupling constants are 
random in a slab of size L'^~^ x A and then periodically 
repeated L/X times along a chosen direction. The distri- 
butions used for the Jij 's vary here from case to case but 
are always chosen so that the interfaces are rough even 
for small lattices sizes, and even in the {100} orientation 
cubic systems. In Fig. |l|, we illustrate the way in which 
the periodic disorder is implemented for the {10} and 
1 11 d irections of a square lattice. As is now well known 
~ll^ , the ground state interface of the system (^) can 
be found exactly using the maximum flow algorithm. We 
have a custom implementation of the push-relabel algo- 
rithm for this problem and using it we are able to find 
the exact ground state interface in Ising systems of size 
one millions sites in about one minute of CPU time on a 
high end workstation. 



III. SCALING THEORY AND NUMERICAL 
RESULTS 

Consider the ground-state interface of a square lattice 
in which the bond disorder has period two in the {11} ori- 
entation (e.g.. Fig. |l|(a)). It is obvious that the interface 
is highly degenerate, as the ground state interface may 
start in any of L/2 equivalent positions. Consider now 
starting to create a ground state interface from the left 
side of Fig. 0(a). To minimize the interface energy one 
chooses the weakest bond. Having chosen this weakest 
bond, the interface crosses this weakest bond and chooses 
the weakest bond in the next column. This process of 
choosing the weakest bond continues across the sample 
and, for period two, the random walk so generated gives 
the exact ground state. The reason this ground state 
is exact is that at each step, all of the possible random 
bonds in each column are tested (there are only two!). 
Thus in this limit, w ~ L^/"^ as for a random walk. In 
contrast, if the period diverges, the model returns to the 
random bond Ising universality class (or equivalcntly the 
directed polymer (DP) in a random medium) for which 
For finite A, we expect that the interface will 
seek to optimize its global wandering until the roughness 
reaches the wavelength of the periodicity 1 19|. After that 
it has exhausted all possibilities and then returns to a 
random walk behavior. We thus have. 



w(L,A) 



w < A, 
w > A. 



(4) 



A natural scaling form based on these limiting behaviors 
is, 



(5) 



where the scaling function /(z) for the roughness has the 
asymptotic behavior. 



const 

^-1/6 



z < 1, 

2 > 1. 



(6) 



Tests of the asymptotic behaviors (^) and the scaling 
function (||) and the results are presented in Figs. || and 
^ for the {10} orientation. It is seen that the predictions 
of the scaling theory are nicely confirmed. Similar results 
were found for the {11} orientation, too. 

We turn now to the behavior of random surfaces in 
(2 -|- 1) dimensions. There, renormalization group (RG) 
techniques have been applied to the random-phase sine- 
Gordon model |p|-|rT|, to random bond interfaces and to 
fairly general models of periodic elastic media. Numeri- 
cally, exact maximum-flow-minimum-cut and minimum- 
cost-matching algorithms and Monte Carlo methods 
pof) have been used. In the random substrate problem, 
there is a low temperature "super-rough" phase where 
~ ln^{L), while in the random manifold problem, the 
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surface roughness is found to behave as w ~ L^^^ , where 
Crb = 0.42 ± 0.01. The quahtative reasoning expressed 
in the first paragraph of this section also applies to higher 
dimensions, so that we expect the behavior of Tip to be in 
the random substrate universality classes at long length 
scales w > \, while the random manifold universality 
class is dominant at short length scales w < X. The lim- 
iting behaviors in dimension (2+1) are then, 

We thus expect, 

-(i,A)^LC«v(^), (8) 

and that the scaling function in (2 + l)-dimensions is 

, , J const z < 1, , . 

^^''> ^ \\uz/z^''^ Z>1, 

with the scaling parameter z — L j I ^"-'^ . The asymp- 
totic behaviors of Eq. (|^) are illustrated in Figs. |^(a) 
and (b) for interfaces in the {100} orientation. The loga- 
rithmic asymptotic behavior is clearly confirmed in Fig. 
§(a), but the random manifold behavior is still strongly 
effected by finite size effects. This is understandable as 
large system sizes are necessary to see the asymptotic 
random manifold behavior, even in the A — >• oo limit 
Though finite size effects are clearly evident in the 
scaling plot of Fig. ^c), the data collapse at large A is 
quite satisfying. It is clear that the random substrate 
(Bragg glass) universality class [^,^ is dominant at 
large enough length scales. We have tested the behavior 
in the {111} orientations and find that {111} interfaces 
behave in a similar manner. 



IV. CONCLUSIONS 
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We have studied the scaling behavior of an elastic man- 
ifold in the presence of a periodically repeated "strong" 
bond disorder. We find that in (1 + 1)- and in (2 + 1)- 
dimensions, and at long distances, the periodicity is rele- 
vant so these interfaces are in the random substrate uni- 
versality class. This is to be contrasted with an interface 
in a system with a periodic potential and with random 
disorder. In the latter problem the periodic potential is 
claimed to be irrelevant on long length scales in (1 + 1)- 
and (2 + l)-dimcnsions for any disorder [0-H, though at 
weak disorder numerical work on {100} orientation cubic 
lattices indicate a strong tendency to order due to lattice 
effects |,|l|. 
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FIG. 1. An example of interface in a random substrate 
problem, with period A = 2: (a) in {11} orientation; (b) the 
{10} orientation. The dotted lines, • • •, describes the lower 
energy bond of the two bonds (in the system of period 2), 

while the dashed line, , describes the higher energy bond. 

A minimum energy path through each system is indicated 
with a thick solid line. 
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FIG. 2. (a) The roughness (w) of manifolds divided by 
the wavelength of the periodicity (A) vs. normalized sys- 
tem size {L/X^^''), where C = (dp = 2/3, for {10} oriented 
(1 + l)-dimensional systems. The random bonds are from a 
uniform distribution with strength AJij^±/Jo = 1 in the per- 
pendicular (z) direction, and AJij_||/Joj| =0.1 in the parallel 
(x) direction in all layers in order to break the degeneracy. 
>^o,||/^o ~ 0.2. The number of realizations = 200 for each 
wavelength A £ [10, .., 160] and system size e [20^ - 1280^]. 



The solid line, — , has a slope ( — (dp 
line, , has a slope ( — (rw ~ 1/2. 
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FIG. 3. The scaling function f(z) — w/L'' of the rough- 
ness w{L,\) vs. scaling parameter z = L/\^^^, where 
C = (dp = 2/3 for the same data as in Fig. |2| The solid 
line has a slope of (rw — (dp = —1/6. 
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FIG. 4. (a) and (b) The data-collapse, w/X vs. L/A^^^ 
C = Crb = 0.42, for the roughness of (2 + l)-dimensional 
{100} oriented systems. The random bonds are from uni- 
form distribution with AJij^±/Jo — 1 in the perpendicular 
[z) direction and constant Jijj|/Jo = 0.2 in the parallel {x,y) 
direction. The number of realizations A'^ = 200 for each wave- 
length A € [3, ..,13] and system size € [6^ - 90*]. (c) 
The scaling function f{z) = w/L'^ of the roughness w{L,X) 
vs. scaling parameter z = L/A^^'-. Finite-size effects with 
logarithmic corrections are visible as a curvature for small L. 
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